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Recent experiments claiming formation of quantum superposition states in near macroscopic sys- 
tems raise the question of how the sizes of general quantum superposition states in an interacting 
system are to be quantified. We propose here a measure of size for such superposition states that 
is based on what measurements can be performed to probe and distinguish the different branches 
of the state. The measure allows comparison of the effective size for superposition states in very 
different physical systems. It can be applied to a very general class of superposition states and 
reproduces known results for near-ideal cases. Comparison with a prior measure based on analy- 
sis of coherence between branches indicates that significantly smaller effective superposition sizes 
result from our measurement-based measure. Application to a system of interacting bosons in a 
double-well trapping potential shows that the effective superposition size is strongly dependent on 
the relative magnitude of the barrier height and interparticle interaction. 



(N 

> 



\0 
O 

Oh 



I. INTRODUCTION 



Despite quantum mechanics being one of the most 
sweepingly successful theoretical frameworks in the his- 
tory of physics, there has always been and still appears 
to be a great deal of unease and confusion about some of 
its fundamental concepts and consequences. Most strik- 
ingly, quantum mechanics requires that if the outcomes 
of certain experiments are known with certainty, then it 
will not be possible to predict the outcome of other, in- 
compatible experiments. Instead, the system must exist 
in an indeterminate state, allowing for the possibility of 
several different outcomes of these experiments. In many 
interpretations, this is viewed as the system simultane- 
ously existing in a "superposition" of all the different out- 
comes at once, until an experiment is actually performed 
and an outcome determined. 

This seemingly ghostly state of affairs is perhaps not 
very unnerving in the context of atoms and microscopic 
systems. But, as Schrodinger pointed out in 1935 a 
microscopic system coupled to a macroscopic one would 
inevitably lead to a situation in which even a macroscopic 
living being — in his example a cat — could conceivably 
end up in a state of being neither alive nor dead, until 
an observer actually looks and determines its fate. One 
"solution" proposed by some people uncomfortable with 
this situation, is that there may be some intrinsic "size 
limit" for quantum mechanics, which somehow prohibits 
nature from putting macroscopic systems into this kind 
of counter-intuitive superposition (see e.g., Ref. 0] for 
a review). Although one may doubt such a proposal or 
question the need for it, it does deserve to be investigated 
whether it can be formulated in a precise enough way 
to be tested experimentally, especially given claims in 
recent years that "Schrodinger cat" states have been or 
can be produced in more or less macroscopic systems 



In order to investigate any possible size limits to quan- 
tum mechanics experimentally, one must of course have a 
reasonably clear definition of what the size of a system in- 
volved in quantum coherent behaviour is. In this paper 
we will investigate systems described by cat-like states 
that can be generically written as \'^) oc \A) + \B), where 
1^) and \B) are macroscopic or mesoscopic states that 
are distinguishable to some extent. The task is to define 
a measure of how "large" this quantum superposition is 
in terms of the constituent subsystems. Each of these 
notions will be made more precise in the course of this 
paper. We explicitly seek a measure that is independent 
of the physical nature of the subsystems and that can 
therefore be used to compare the effective size of cat-like 
states realized in very different physical situations, e.g., 
Bose Einstein condensates and superconducting current 
loops. 

This question, which could be succinctly phrased as 
"how big is Schrodinger 's cat" for a given system in a 
particular quantum superposition state, has been asked 
in several earlier papers [3^^. By size we mean the num- 
ber of effective independent subsystems that can describe 
the superposition (we will discuss in more detail what we 
mean by these notions in Section |TT| . One "ideal" N- 
particle cat state, for which the answer would be N, is 

a GHZ-state of the form |*) = 2-^/^ f |0) 
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where |0) and |1) are any pair of orthogonal one-particle 
states. Hardly any states realizable in the laboratory are 
of this idealized form however, and we therefore seek a 
measure that can quantify the size of more general states 
that are still recognizable as generic cat-like states but 
that may be very different from the ideal form. The par- 
ticular case of a generalized GHZ-like A^-particle state of 

the form |*> = R-^ (jcl^of'^ + , where |(^o) and 

{(pi) are non-orthogonal one-particle states, was studied 
in 8] with two independent approaches, one based on 
the stability with respect to decoherence, and the other 
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on the amount of distillable entanglement. In the hmit 
of highly overlapping states, where \{(/)o \ = 1 — 

for some e <C 1, Diir et al. found that both decoherence 
and distillable entanglement measures of the "effective" 
number of degrees of freedom participating in the super- 
position, n, yielded an effective cat size of n ~ Ne^ 18|. 
These two measures were specific to the form of the non- 
orthogonal GHZ-like states and it is not obvious how to 
apply them to arbitrary superposition states. Another 
motivation of this paper is thus to derive a measure of 
the effective cat size that can be applied to superposi- 
tions \A) + \B) having completely general forms of the 
states \A) and \B). 

The rest of this paper is divided into four parts. In 
section 2, we present a measure of effective "cat size" 
for general binary superposition states that is based on 
the notion that the "cattiness" of a superposition state 
should depend primarily on how distinguishable the two 
branches of the state are. This measure is based fun- 
damentally on measurements, and thereby differs from 
earlier measures that have tended to be based on the 
mathematical form of the state. The new measure is 
thus potentially more useful for experimental implemen- 
tations. In section 3, we apply the measure to a system of 
bosons in a two-mode description. In section 4 we con- 
nect the results from section 2 with realistic numerical 
Monte Carlo simulations of Bosons with attractive in- 
teractions trapped in a double-well potential. Section 5 
summarizes and indicates future directions of research. 



II. IDEAL CATS AND EFFECTIVE CAT SIZES 

In this Section we will give a definition of the size of a 
cat-like state of an object, I^E*) = \A) + \B). We will con- 
sider that the object is formed by N subsystems, and our 
measure of effective size will then range between and N, 
analogously as in Ref. fE]. However, in contrast to that 
work, the quantity we introduce here will measure how 
(macroscopically) distinguishable the states |^) and \B) 
are. The main idea that we want to capture with this def- 
inition is the following: how many fundamental subsys- 
tems of the object do we have to measure in order to col- 
lapse the entire state into a single branch corresponding 
to one of the two states \A) or \B), and how many times 
larger than this number is the entire system? By "fun- 
damental subsystem" , we mean something that in some 
sense can be taken as a fundamental building block of 
our system, e.g., single particles or something similar. It 
is by no means always clear what one should consider the 
fundamental building blocks of a given physical system 
(molecules, atoms. Cooper pairs, electrons, quarks...), 
and we will not attempt to make a definitive definition of 
what such building blocks should be, if this is even pos- 
sible. However, our measure will be based on how many 
measurements must be carried out to perform a specific 
task, namely to collapse the superposition state into one 
branch or the other. A reasonable qualitative definition 



would therefore be that a fundamental subsystem is the 
smallest subsystem that one could in principle measure 
in some experimental context and which would provide 
information that could help distinguish one branch from 
the other. For a BEC experiment one could in principle, 
e.g., scatter light from single atoms, making single atoms 
reasonable candidates for fundamental subsystems. Our 
measure thus will depend on the experimental situation 
and the relevant size and energy scale, something which 
probably must be expected if one wishes a measure that 
does not involve Planck-scale physics. For the remainder 
of this paper, even though relevant fundamental subsys- 
tems may not always be something that can reasonably 
be called particles, we will use the terms "particle" and 
"fundamental subsystem" interchangeably, and this con- 
cept plays an important role in our measure. More specif- 
ically, the question we ask to define our measure is the 
following: What is the maximal number of disjoint sub- 
sets that one can constitute from the N particles such 
that by measuring all particles in any given subset one 
can cause the superposition state to collapse into one 
of the branches \A) or \B) to a high degree. A mea- 
surement that causes such a collapse is equivalent to a 
measurement that with high probability lets us deter- 
mine correctly whether a system is in state \A) or \B) if 
we are given a system which is definitely in cither one of 
these two states, but we do not know which one. We em- 
phasize that the latter situation is clearly very different 
from having a system which is actually in a superposition 
\A) + \B). But since a measurement which collapses the 
superposition state is idential to one which is capable of 
distinguishing between the two branch states (assuming 
an ideal measurement with no classical noise), we shall 
often use the latter picture in the discussion below 

It is not difficult to write a mathematical definition 
which expresses our measure as formulated above. How- 
ever, in practice it may be quite difficult to calculate this 
for general superpositions, since for a given accuracy one 
has to optimize the number of subsets over all possible 
partitionings of the N particles. Thus, we will use an 
alternative definition that also captures the above con- 
cepts but is simpler to evaluate, particularly for states 
possessing permutation invariance. 

Definition of cat size. Given an object composed of N 
subsystems and < S ^ 1, we define the cat size of a 
state l^-) oc \A) + \B) with \\ \A) \\ = || \B) \\ = 1, to a 
precision S, by 

C5(*) :-iV/n„,in, (1) 

where nmin is the minimum number of particles one has 
to measure, on average, in order to distinguish the states 
\A) and \B) with probability greater than or equal to 1-5. 

In order to determine Cs{'^) we can proceed as fol- 
lows. We begin with 1-particle measurements {n — 1). 
For each particle k we calculate the optimal probability of 
being able to distinguish \A) and \B) by measuring just 
this particle and average this probability with respect 
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to fc. If the resulting average probability is larger than 
1 — S, then nmin = 1 and hence Cs{'i') = N. If not, we 
then go on to consider all possible sets of two particles, 
{j,k), determining the corresponding optimal probabil- 
ity of distinguishing \A) and \B) by measuring these two 
particles. If after averaging this probability with respect 
to J, k we obtain an average probability larger than 1 — S, 
we have n^in = 2 and hence Cs{'^) — N/2. If not, we 
repeat the procedure with measurements of an increas- 
ing number of particles until we reach a value of n^i^ 
for which the averaged probability of successfully distin- 
guishing the two branches is for the first time larger than 
1 — S. If this happens only when all particles are mea- 
sured, then Timin = N, and the cat size is Cs — 1. If 
even measuring all N particles still fails to distinguish 
the two branches to the desired precision 1 — 6 , then 
n,nin and hence the cat size Cg are essentially undefined. 
For simplicity, we will define the cat size to be zero in 
this situation. 

Thus, the only ingredient we need in order to de- 
termine the cat size is the maximal probability to be 
able to distinguish two states |^) and \B) by measur- 
ing only a given subset of the total system (or using 
some similarly restricted set of measurements, as we will 
see in Section HlH). We now briefly discuss this proba- 
bility. For more thorough and general discussions, see 
Refs. [1, [13, Using a generalized quantum mea- 

surement, i.e., a POVM (positive operator valued mea- 
sure) in which the outcome described by POVM 
element Ea is taken to indicate that the system is in 
state \A) and the outcome Eb that it is in \B), then 
given equal prior probabilities for each state (i.e. equal 
weight for the two branches of the superposition), the 
probability P of inferring the correct state from a single 
measurement is 



P 



1 



tr {pAEA)+tr (pbEb) 



(2) 



where pA — \A) {A\ and ps = \B) {B\ are the density ma- 
trices of the two states. If we now restrict ourselves to 
measure only a subset of n particles, then the measure- 
ment outcomes are given by POVM elements E^\ i?^"' 
that act non-trivially only on these n particles, acting 
as the identity on the remaining N — n particles. The 
probability of successfully inferring the state is then 



P = - 



tr 



tr 



PA E^ 
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(N-n) 



tr pB E 
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with p^'' = 



trj^-nPA and p^-* = tvN-nPB the corre- 
sponding n-particle reduced density matrices (n-RDMs) . 
(trjv-ri denotes the trace over all particles except the n 
particles being measured.) The maximum probability for 

(n) 

successfully distinguishing two density matrices Pa and 
p^'' will then be given by an optimal POVM, which is 



known to be a projective measurement in the eigenbasis 
of the operator p^"^ - [3, [H 



(4) 



Here ||X|| = tr|X| is the trace norm, i.e., |Ai|, with 
Ai the eigenvalues of X. 

Several remarks are in order here. 

(i) We have based our working definition here on the 
average probability over all equal size subsets being 
larger than 1 — 5. One could alternatively have em- 
ployed a requirement that the minimal probability 
is larger than 1-5. Also, as mentioned above, at 
the cost of introducing a great deal more computa- 
tional expense, one could replace the average over 
equal size subsets by the optimum partition over 
all possible subsets. 

(ii) Although it should be clear from the notation, we 
note that, as defined, our measure applies only to 
pure quantum states, not mixed states. Defining 
a cat size measure for mixed states is complicated 
by the fact that there is no unique way to decom- 
pose a mixed state density matrix into a convex 
sum of pure states, so that, e.g., a mixture of com- 
pletely separable states could also be written for- 
mally as a mixture of very cat-like states. Any 
cat size measure applicable to mixed states would 
therefore have to weight the cat size quite heavily 
with the purity of the state. We will not pursue 
such an extension of our measure in this paper. 

(iii) For states that are symmetric with respect to per- 
mutations, for a given number of measured parti- 
cles n it suffices to consider only a single subset, 
since all subsets give rise to the same probability 
because of symmetry. This results in a consider- 
able gain for computational studies with large N 
and will be analyzed in detail for bosonic systems 
in the remainder of this paper. 

(iv) We have assumed that we can perform collective 
measurements on a subset of n particles. However, 
we can also consider the situation in which only 
individual single-particle measurements are per- 
formed. In some cases the calculation could then 
be highly simplified, since we would have to con- 
sider only single-particle reduced density operators. 
This situation appears well suited to bosonic sys- 
tems and will be analyzed further in Section IIII Bl 

(v) Given a state in which \A) and \B) are not 
specified, there are many ways of selecting the two 
branches, and these may give rise to different values 
of the measure. Thus, when we talk about the size 
of a cat state, we must always specify what are the 
branches A and B. Furthermore, application of the 
measurement-based cat size defined above requires 
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that the two branches have the same norm. If the 
norm of the two branches are different, i.e., 



1*) cx \A)+g\B) 



(5) 



with < g < 1, we expect that Cs{'^) must be 
multiphed by a factor that interpolates smoothly 
between a value of zero when g = Q and a value of 
unity when \g\ = 1. This factor can be determined 
by recognizing that the general superposition for 
general g can always be distilled to the equal super- 
position \g\ — 1 by generalized measurements Q, 
yielding an effective cat size that is reduced by the 
associated probability. For the state of Eq. ([5|) , one 
can perform a measurement using the operator 



^1 



g\A){B^\ , \B){A, 



{B± I A) {A^ I B) 



(6) 



and complement this with any other measurement 
operator A2 such that Ei = A\Ai and E2 = ^2^2 
form a POVM, i.e., Ei + E2 ^ 1. and 
are any states that are orthogonal to \A) and \B) re- 
spectively. If one obtains the outcome correspond- 
ing to Ai , then the state after the measurement will 
be the equal superposition state l^*) cx |A) -I- jS). 
The probability for this to happen is 



Pa 



{2+{A\B) + {B\A))\g\^ 
l + {A\B)g + {B\A) g* + \g\' 



(7) 



Thus if the norm of the two branches are different, 
we can take the effective cat size to be pgCsi'if), 
where g is the smaller of the two norms. Note that 
if at least one particle separates out in each of the 
branches \A) and \B), i.e., if j^*) can be written 
in the form \a) \A^~^') + \b) for some one- 

particle states I a), \b) and {N — l)-particle states 
then the distillation can be ac- 
complished using only a local single-particle mea- 
surement, namely 



Ai 



9W){b\ , \b){a^\ 



{b 



±1 a; 



(a± I b) 



(8) 



The probability of obtaining the outcome Ai is the 
same as m Eq. dZl), with \A) replaced by \a) | A^-^). 

(vi) In order to calculate Cs{'if), we can calculate P = 
1 — Pe , where Pe is the probability of error in dis- 
tinguishing the two states, and then find the value 
of n for which P > 1 — 6. In Section IIIII we will 
show plots of Pe rather than P, since these bet- 
ter illustrate the scaling of the error with n. For 
large N values, in some situations we can also solve 
P — 1 — (5 to obtain a continuous value of n (see 
Section UnBl). 

(vii) Our approach of asking how many subsystems a 
system can be divided into such that each one alone 



suffices to distinguish the branches of a state, has 
some similarities with the concept of redundancy, 
introduced in a different context in . There, the 
redundancy of a piece of information about a quan- 
tum system is defined as the number of fragments 
(partitions, in our terminology) into which the envi- 
ronment can be divided such that this information 
is contained in every one of the fragments. This is 
used in [l^ to probe how objective a certain piece 
of information about a quantum system is, since 
information that has a high degree of redundancy 
can be obtained by many observers independently 
through measuring different parts of the environ- 
ment, without disturbing the system itself or each 
other's measurements. 

(viii) Finally, we note that our measure does not look at 
the physical properties of the object, such as mass 
or spatial dimensions, but rather at the number of 
components. For example, with this measure a very 
massive elementary particle can have a cat size of 
1 at most. 

We now give some examples of the cat size for sim- 
ple superposition states, calculated using the above for- 
malism in a two-state basis. Suppose we have a sys- 
tem consisting of a macroscopic number N of spin- 1/2 
particles. First, consider the ideal GHZ states |^'±) := 
(|0)'*^ ± jl)*^^). Here only one particle need be mea- 
sured to distinguish the two branches with certainty, 
i.e., the one-particle reduced density matrix (1-RDM) al- 
ready gives P = 1, and hence rimin = 1 and Cs = N 
for all S. Now consider the linear superposition state 
|*> = ^ (|1'+> + |*_)) = |0)®^. This is also a super- 
position of two distinguishable (orthogonal) macroscopic 
quantum states, but here all N particles must be mea- 
sured in order to distinguish the two branches. The n- 
RDMs for |*+) and are identical for all n < N, 

so P = unless n = iV, in which case P = 1. Hence 
'^min = N and the cat size is equal to 1, as expected since 
the state is equivalent to a product state. As a final ex- 
ample, we apply our measure to the non-ideal state with 
non-orthogonal branches that was studied in Ref. [1], 
namely |*) {\0f^ + \^f^) with |(0 | </>) p = 1 - e^, 
where e ^ 1. Here, the two branches |0)®^ and 
are separable states, and their respective n-RDMs are 
therefore equal to density rnatrices of pure n-particle 

respectively. In general. 
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states, namely |0) and 
for any quantum system and any pair of states |a) and 
1 5) with I (a I 6) p = c^, we can write the correspond- 
ing density matrices in a two-state partial basis defined 
by I a) and \aj_), where \aj_) is the state orthogonal toja) 
but contained in the subspace spanned by \a) and \b) [37| . 
Specifically, writing \b) — c\a)+s \a±) with |cp-f |sp = 1, 
we have 



Pa~ Pb = 



-sc 



-sc 

lo|2 



(9) 
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Using Eq. (|4|), we find that \a) and \b) can be successfully 
distinguished with probability P=i(l + |s|). Defining 



\a) = |0)^" and \b) ^ 
mum success probability 



we then obtain the maxi- 



P = ^(l + v/i-(i-^')") 



(10) 



for distinguishing jO)**^ and 



\0N 



using ri-particle mea- 



surements. Requiring this to be greater than 1 — 6, where 
S is the desired precision, results in a value of nmin given 
by 



macroscopic cat states, are performed on systems of iden- 
tical particles. These include photon states [ij], super- 
conducting current loops [1, 01, spin-polarized atomic en- 
sembles 5] and Bose Einstein Condensates Ts']. Cat 
states of bosonic particles allow some simplification of 
the proposed measure of effective cat size, since making 
use of the permutation symmetry reduces the size and 
number of the n-RDMs to be analyzed. We consider here 
a generic form of cat state wavefunction that generalizes 
the ideal GHZ state 



log(4J-4J^) 
log(l - e^) 



(11) 



where [. . .] denotes the ceiling function, i.e., the nearest 
integer above the value of the argument. For e and S 
small, this results in Cs — A^/?T-min — Ne^ /{— log(J)), in 
agreement with the Ne^ scaling found for these states in 
Ref. H. 



\GHZ_ 



N 



V2 



(12) 



III. CAT STATES IN BOSONIC SYSTEMS 

Most experiments involving quantum coherence in 
more or less macroscopic systems, including potentially 



to situations described by a superposition of non-ideal 
GHZ-like states in which the single particle states are 
non-orthogonal. In particular, we consider states of the 
form uM 



1^-) cx / def{e) (cos 6* + sin 61 6^) 



N 



eb^f |o) 



smWa' + cos 



(13) 



where the operators and 6^ create two orthogonal 
single-particle states. For fixed values of 0, the inte- 
grand of Eq. corresponds to ground states of a two- 
state BEC with attractive interactions, found in [l3| us- 
ing a two-mode approximation and an extended mean- 
field calculation. In this section we will illustrate the 
effects of f{9) for various values of its mean and vari- 
ance. The two branches of the superposition '^'[^'''^ 

''■^''^ are thus defined here by a superposition 



and 



B 



of states 



^(fHfl)) and 



{9) ) that are themselves 



non-ideal GHZ-like states of variable orthogonality de- 
fined by the angle 9. In the notation above, 9 = and 
tt/2 correspond to perfect orthogonality of the single- 



particle states 

(sin 9 



(cos 61 flt + sm0b^) |0) and 



-cos 



>b^)\0) (with 



|0) 



and 



(t)^g\0)'^ = &^ |0), switched for 6* = 7r/2), 6* = 7r/4 



corresponds to complete overlap (with both 
(7r/4)\ equal to 2-^/^{a^ + b^ |0>), 



6^(^/4) 



and 



= -7r/4 

also corresponds to complete overlap but with differ- 



ing overall sign ( 
2-i/2(at _ 5t) |o)) 



(-V4)^ 
7r/2 corresponds to or- 



'(-7^/4), 
and : 

thogonality again but with a factor of —1 for each of 
the states relative to 9 = 7r/2. The extent to which the 
two branches can be delineated is clearly dependent on 
the amplitude function f{9) that controls the amount 
of spreading of each branch. The form of the spreading 
function f{0) will depend on the details of the physi- 
cal realization of the macroscopic superposition, as will 
the values of the angle 9. This generalized juperpo- 
sition reduces to the form employed in Ref. 17^ when 
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f{9) = S{9 — 9o) for some Oq dependent on the parame- 
ters of the Hamiltonian used there, and is in agreement 
with general expectations for the form of macroscopic 
superposition wave functions for superconductors ^ . In 
Section |V] we analyze the form of f{6) appropriate to a 
cat state formed from a BEC trapped in an external dou- 
ble well potential. Numerical calculations for attractive 
Bose gases have shown that the competing effects of tun- 
neling between modes and interactions between particles 
can be taken into account by letting f{9) be a Gaussian, 
the shape of which is determined by the ratio of tun- 
nelling and interaction energies Note that while 
Eq. ([13]) is implicitly a two-mode wave function, this form 
can readily be generalized to multi-mode superpositions. 



A. Calculation of effective cat sizes for 
superpositions of non-ideal states 

We can first give some qualitative expectations for the 
effective size of this superposition state when TV becomes 
large. There are two factors that will reduce the effective 
size below that of the ideal GHZ state, N. Firstly, for 
values oi 9 ^ 0, the two branches of |^) in Eq. are 
not orthogonal, and hence not completely distinguish- 
able. As shown explicitly in Section [IT] above, our mea- 
sure therefore gives a cat-size for this state that is smaller 
than iV, in agreement with the results derived previously 
in Ref. Second, if the amplitude function f{9) de- 
viates from a 5-function, the inner product between the 
two branches will not approach zero even in the limit 
N —^ oo. Hence there will always be a finite minimal 
probability that we will not be able to distinguish the 
two branches, even in the thermodynamic limit and even 
if all N particles are measured. Eventually, if this irre- 
ducible overlap between the branches is large enough, the 
division into two different branches becomes meaningless. 
The effect of this second factor has not been investigated 
before, but is essential to investigate for understanding 
macroscopic superpositions in realistic physical systems. 

To make quantitative calculations for states of the form 
of , it is convenient to first make a change of basis as 
follows, 

c=^{a + ib) d^^{b + ia) (14) 

ct = i=(at-^6t) dt = -L(6t_^at), (15) 

so that the integrand components of the two branches in 
Eq. (fT3l) become 



1 



(0) 



1 



(e*''(i^+ze-*^c^)^|0). (16) 



When measuring indistinguishable bosons, we ob- 
viously cannot pick out n specific particles to make 



an n-particle measurement as described in the dis- 
cussion in Section [Hi For indistinguishable particles, 
the Kraus operators [lO] describing the effect of any 
measurement outcome have the form, e.g., ^["^ = 
c).^*^ "*"aiiai2 • • • , where i denotes a single- 
particle state, with corresponding POVM elements 



•in t t 
a'- • • • a'- a 



{i},{3} 



(17) 



Here k labels the outcome and the superscript (n) spec- 
ifies the number of particles on which the operator acts 

Eq. (fT7| gives us the probability 



M. = tr I*) (*| E, 



(») 



tr 



(18) 



for a given outcome when the system is in state 

1 4'). Here is the matrix given by the coefficients 



1^ , X m/{N ~ny., and 



(„)yi»2-»„ ^ {N-ny. 



X (*|aL •••4al%-iaj2 •••aj„ l*> (19) 

is the n-particle reduced density matrix, or n-RDM, of 
the bosonic system in second quantized form. The combi- 
natorial factors here are introduced so that p*^"-' will have 
trace 1. Furthermore, since p*^"-' is symmetric in both all 
upper and all lower indices, we can index p*-"^ by mode 
occupation numbers k and I. The resulting symmetrized 
matrix acts on a vector space which is equal to the full 
vector space projected onto a symmetric subspace |21j]. 
Denoting the symmetrized RDM by p^"^ we obtain: 




(20) 



where the index k refers to the number of creation op- 
erators equal to and / to the number of annihilation 
operators equal to c 39]. With this definition, the sym- 
metrized n-RDM p^"' has the same nonzero eigenvalues 
as p^") and can therefore be used in place of p^"^ for the 
calculation of effective cat sizes. 

This projection onto the symmetric subspace results 
in a significant reduction in dimensionality, permitting 
calculations to be made for values of n up to several hun- 
dred. Matrix elements of pI^*'' —p^"* are readily calculated 
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for general forms of the amplitude spreading function 
f{9) (see Appendix E]). A key component of these matrix 
elements are inner products between the states \4>a b(^)) 
at different values of 0, which yield factors of cos^{0 — 6') 
and sin^(6' + 9'). For large values of N these functions 
can be approximated by delta functions. This simpli- 
fies the resulting integrals but removes any explicit A^- 
dependence from the result (see Appendix]^. The ma- 
trix pI^' — p'"g ^ is then diagonalized and Eq. ^ evaluated 
to obtain the maximal probability of successfully distin- 



guishing 



*f and 



with an n-particle measure- 



ment. The effective cat size Cs is then obtained by de- 
termining the minimum value of n such that P > 1 — 5, 
according to Eq. ([T]). When using the delta function 
approximation for large N, since the total number of 
particles is unspecified, we evaluate the relative cat size, 
Cs/N = l/n,„in. 

Figures [1] and [2] show the results of calculations for a 
gaussian amplitude spreading function 



(O-Oq)-' 



(21) 



This form is convenient for a systematic analysis of the 
behavior of effective cat size with spread and overlap of 
the two branches since all matrix elements are analytic 
(see appendix The range of 6*0 should be from — 7r/2 
to +tt/2 in order to encompass all relative phases and 
degrees of overlap/orthogonality. Superposition states 
characterized by cr = possess zero spread and reduce 
to the non-ideal states studied earlier in Ref. 8] that 
are characterized by the extent of non-orthogonality for 
6*0 > 0. Figure [1] shows the error probability Pe = 1 — 
P, plotted on a logarithmic scale as a function of n, for 
various values of the spread function parameters Oq and 
a. We show Pe rather than P, since the former allows a 
clearer analysis of the differences between results for a = 
and for (T 7^ 0. The relative effective cat size Cs/N = 
^/nrnin resulting from these probabilities is plotted as a 
function of and a for several different values of the 
precision parameter d in Figure [2l 

Figure [T] shows that while for all values of the parame- 
ters ^0 and a there is a generic increase in the probability 
P for distinguishing the two branches of the cat state as 
n increases (i.e., a decrease in the error probability Pe), 
the nature of this decrease is strongly dependent on the 
actual values of 9q and a. For a — 0, the error is due 
entirely to non-orthogonality, as discussed in Ref. {§] and 
Section |TT1 Here, when = the generalized super- 
position reduces to the ideal GHZ state and the error 
probability is zero, independent of n (not shown in the 
bottom right panel since the logarithmic scale cannot ac- 
commodate Pe — 0). When 0, the non-orthogonality 
makes the success probability increase more slowly with 
n, and hence the effective cat sizes in Figure [D become 
smaller as 6*0 approaches the value ±7r/4 at which the 



two branches 



5- 



and 



^ / 



particular, for strong overlap. 



overlap completely. In 



e <C 1 (outer limits of on ct = axis) , we verify that the 
relative cat sizes are in accordance with the asymptotic 
scaling ~ established in Section [III This effect of non- 
orthogonality also acts when cr > 0, with the relative cat 
sizes also dropping off away from = 0- However now 
there is an additional decrease, due to the branches of the 
cat state getting "smeared out" and overlapping more as 
the width parameter a increases. For all cr, we see that 
the effective cat size is largest for 6o = 0, where the two 



branches 



and 



are orthogonal. 



Detailed analysis of the dependence of the error prob- 
ability Pe on the width parameter cr provides additional 
information. When cr = and Oq ^ 0, consistent with the 
scaling shown in Section|II]the error probability decreases 
exponentially with n and asymptotically approaches zero 
as more particles are measured (solid blue lines in top 
right and bottom panels). However, for cr > 0, we see 
that the decrease in the error probability is slower than 
exponential. In fact it appears to never approach zero but 
is instead bounded below by some finite value, implying 
that the success probability is bounded away from unity. 
This derives from an important feature of this gaussian 
amplitude function f{6) that is illustrated by comparing 

the overlap between ^'^''^ and ^'^''^ for different val- 

For example, at cr = 0, 7^ 7r/4, the 



ues of cr and 60 ■ 
inner product between 



5- 



^ / 



and 



^^^\ goes to zero 



as A^ ^ 00, so that the two branches become orthogo- 
nal in the limit of an infinite number of particles, and 
one can therefore always tell them apart with arbitrarily 
high certainty by measuring enough of the particles (solid 
blue line). However, for a > 0, the overlap approaches a 
finite value as A^ — s- 00. In this situation it is not always 
possible to distinguish the two branches within a given 
precision, regardless of how many particles are measured 
- even for n — N. This implies that nmin is undefined for 
these extreme cases. As noted in Section [TTl we formally 
define Cs = in these situations, with the additional un- 
derstanding that l^*) is not really a meaningful cat state 
at all here. 

This behavior for cr > is consistent with the fact 
that the two branches 



0f )(0)) and 



b^j^\0) ) can 



interchanged, either by transforming 6 
< 6 < tt/2, or by first transforming 9 
then changing sign, for > 6 > — Tr/2. 
amplitude spread function f(0) has support both inside 
and outside the region — 7r/4 < 9 < +7r/4, some of 



be 

7r/2 - 6*, for 
> -7r/2- and 
Thus when the 



contributes to both branches 



and 



and 



the state cannot be split into two disjoint branches. Us- 
ing Eqs. (jAl[) -(|A3 |) . it is also easy to see that for cr = 0, 

> when N ^ 00, so that the branches 



(N) 



{N) 



become orthogonal and distinguishable in the thermo- 



(AT) 



(N)\ 



I 



ap- 



dynamic limit, whereas for cr 7^ 0, x^-a 

proaches a finite minimum value. This is the physical 
reason why two strongly overlapping branches cannot be 
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FIG. 1: (Color online.) Error probability, Pe = 1 — P, for distinguishing the two branches of the generalized cat state 
superposition Eq. (|13|l when characterized by a gaussian amplitude spreading function f{6), for various values of the gaussian 
parameters 6q and a. 



distinguished to arbitrary high precision ((5 — > 0), even in 
the hmit N,n —i- oo. Detailed analysis of the support of 
the amplitude spread function will thus be very impor- 
tant for realistic estimates of cat size in physical systems 
involving superpositions of non-orthogonal states. 

This difference in behavior of success probability scal- 
ing for a = and for a > has a large effect on 
the effective cat size. Figure [5] shows the effective rel- 
ative cat size Cg/N for four different precision values, 
5 = 10-2, 10-4, 10-^ and 10-^°. It is evident that if S is 
sufficiently small, the effective cat size does not depend 
too heavily on the exact value of 6 when a — 0. This is 
to be expected, since 1 — P decreases exponentially with 
n when cr = 0, and hence rtmin will only be proportional 
to log{l — P). However, when cr > 0, we see that the 
cat size can be significantly reduced or even vanish for a 
given system as we decrease the desired precision S. This 
illustrates the point made above, namely that states with 



a > become increasingly poor cat states as a increases 
and eventually are not cat states at all. It also provides a 
dramatic illustration of the general fact that the degree 
to which a superposition state can be viewed as a cat 
state is inherently dependent on the precision to which 
the implied measurements are made. 



B. Estimate of effective cat sizes from 
single-particle measurements 

In all of the analysis so far, we have assumed that any 
71-particle measurements can be made to distinguish the 
branches \A) and \B) of a cat state, including collective 
measurements in entangled bases. In practice, this is usu- 
ally not feasible for large values of n. From a practical 
perspective, it would therefore be desirable to have a def- 
inition of cat size which relies not on general n-particle 






FIG. 2: (Color online.) Relative effective cat size Cs/N = l/rimin as a function of the gaussian parameters 60 and a, for 
several values of desired precision 5. All plots have a resolution of 7r/40 in both 9o and a. Numerical calculations were made 
for n < 100, imposing a numerical cutoff of 0.01 on the value of l/umin- 



measurements, but instead only makes use of measure- 
ments that can be put together from n separate 1-particle 
measurements. 

Allowing only those n-particle measurements that can 
be realized as a sequence of 1-particle measurements 
means that we restrict the corresponding POVM ele- 
ments to be of the form 

^=Em£^i:^4f (22) 

where each e\^'' = A-'^''^j4,-'°'' acts on a single particle k 

*fc •'k '■k 1 

only, and where pi are positive numbers subject to the 
constraint that tr£ < 1. (Note that, unlike the situation 
in Sections |TT] and IIII Al the POVM elements here act 
each on only a single particle, and the superscript index 
(fc) in parentheses therefore labels the particle that each 
operator acts on, not the number of particles it acts on.) 



This means that the POVM elements must be separa- 
ble. Furthermore, to ensure that the measurements can 
be realized as a sequence of 1-particle measurements, it 
must be possible to write express the POVM elements 
in such a way that e'^^'^ only depends on E^^^^ for I < k 
but not for I > k. To find the maximum probability P 
of successfully distinguishing the branches \A) and \B) 
of a cat state using such measurements, we would then 
need to maximize Eq. ([3]) with E^^^ and E^^^ subject to 
the above constraints. Unfortunately, we know of no ef- 
ficient way to do this. In particular, deciding whether a 
given POVM is separable as in Eq. (|22p is known to be 
an A^P-hard problem [2^ . 

However, if we restrict ourselves to a very simple case, 
namely to superposition states where each of the two 
branches of the cat state are themselves product states, 
not only is the optimal measurement strategy using a 
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sequence of n one-particle measurements known, but it 
even performs equally well as the optimal general n- 
particle measurement. To show this we adapt the tech- 
niques used in [23| . In that work, one is given n copies of 
a quantum system, all prepared in one of two states \iPa) 
and \iPb) and asked to tell which one (note that [1^ uses 
0,1 rather than A,B). The joint state of all n copies 



is then either \iPa)^"^ or lips)'^"', and the corresponding 



0n 



density matrix is p[ 



(gin — 



for ^ = A,B. 



One assumes prior probabilities qa and qs — 1 — qA 
that the correct state is \iPa) and IV's), respectively. The 
maximum possible probability of guessing the right state 
would in general consist of making an optimally chosen 
collective n-party measurement (i.e., possibly in an en- 
tangled basis) on the n copies. However, it is shown that 
by measuring only a single copy at a time and choosing 
each measurement according to a protocol that effectively 
amounts to Bayesian updating of the priors qA and qs 
based on the outcome of the previous measurement, one 
can obtain a success probability which is equal to the 
maximum one for a general rt-party measurement. 

In our case, we are trying to ascertain whether a single 
system consisting of N subsystems is in a state \'^a) or 
another state Ivl/s), where these states are known to be 
product states with respect to the N subsystems. We 
can therefore write 



(1) 



(2) 



(23) 



where 'y = A,B and 



/(fc) 

■07 



is the state of particle num- 
ber fc, and we assume that we will measure the first n 
particles. This is equivalent to a generalization of [l^ 
to a situation where not all the copies of the system un- 
der study are the same, but where each "copy" k is in 

one of two states V't-'^''^ for 7 = A or _B, and where 7 
is the same for each k, and the task is to determine the 
value of 7, by only measuring n of the "copies" . We will 
now show that the conclusion of ^] still holds in this 
case, namely that the performing a sequence of n opti- 
mal one-particle measurements with Bayesian updating 
between each measurement gives the same probability 
of success as the best collective n-particle measurement. 
We will use a slightly different approach than us- 
ing 1-particle reduced density matrices instead of single- 
particle state vectors, since this approach is more readily 
generalizable to indistinguishable particles. 

Following the notation of [2^, we will here write the 



states 



and 



of particle k in the branches 



I^^a) and I^b) respectively as 



h \xk) + (-l)"sin6'fc \yk) 



(24) 



where a — for 7 = A and a = I for 7 = S, and 
\x) and \y) are two basis vectors in the state space of 
particle k chosen such that this relation is valid (this 
is always possible). The corresponding reduced density 



matrix with respect to particle k in the {|a;),|j/)} basis 
are then 



Ik) ^ { C0S2 dk 

' \{-lY COS Ok Sin Ok 

_ ( cos^Ok ^^sin26'fe 
-l(zirsin2^fe sin^^fc 



(— 1)" cos^fc sin6'fc^ 
sin^ Ok 



(25) 



(note here that the superscript k again refers to the par- 
ticle to which the RDM belongs, not the number of par- 
ticles described by the RDM, which in this case is just 1.) 
If we now let the probability, prior to measuring particle 



fc, of the state being 



■07 



be gV > 



then the measurement 



which produces the highest probability of successfully 
identifying the correct state, is a projective measurement 
in the basis in which the matrix T^^'> = q\^^ p\^^ — q'l^p'i^ 
is diagonal (0,1131). The conclusion 7 = A is associated 
with the eigenspaces with positive eigenvalues of T^^\ 
while -f = B corresponds to the eigenspaces with nega- 
tive eigenvalues. In the basis {\xk) ,\yk}}, the matrix F''^^ 



p(fc) 



q^ ) cos I 



hgf))sin2A 

(fc)^ ■ 2a I ^ 
9b )sm 9 J 



and is diagonalized by 



U(0fe) 



cos <pk sm (pk 
- sin (pk cos (j)k 



(27) 



with 



(fe) i (fc) 1 
sm2(j)k = '^^ ^'^^ sin26>fc ^ — sin26lfc (28) 

Rk Rk 

(k) _ (k) 

COS 2(j)k = „ cos 2ek 

Rk 



(29) 



Rk^i [qf + q's'f MTIb' COs2 2ek 



{k)Jk) , 



= ^ 1 - 4q^^\^^'^ cos2 26 k 



resulting in eigenvalues 



Ak) 



(30) 



(31) 



The outcome E\ is associated with the eigenspace of 
p(fe) corresponding to the eigenvalue which is the 

first eigenvector in the diagonal basis. In the basis used 
in Eq. (pS)) . we then have 



E^^Uic^kY (J o)c/(0.) 



^ f cos^ (t>k \ sin 2(j)k 



^ sin 20J. sin^ 



(32) 



Combining this with Eq. (P5|) gives us the conditional 
probabilities P{e''^^\^) = tr (^Ej^^pij'^^ of obtaining the 
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(k) 

outcome E)^ when measuring particle fc, given that the 
initial state of the joint system was l^'^): 



of the n particles to be measured, and hence the overall 
probability of success, is equal to 



P(4'=) I A) = 1 + -L (l - 2q'^^^ cos^ 29, ) (33) 



P{E^^'^ IB) = 1 - -L (l - 2q^^^ cos' 20, ) . (34) 



The corresponding probabilities of obtaining Eg '' ~ 1 



-E^'^' are then 



p(4'=)|A) 

PiE^lB) 



1 

2i4 



2q^^^ cos^ 29, 
2q^^^ cos^ 29, 



(35) 
(36) 



Using Eqs. ((33)) and p6)) . the probability of successfully 



identifying the state 



after measuring particle k 



(conditional upon ealier measurements yielding the priors 

q^A^ and q^*'^) is 



P, EE qfp{E'^^>\A) + q'^>P{E''^'\B 
^2 + 2^" 



(37) 



To find the overall success probability of the proce- 
dure, we need to evaluate what the posterior probabil- 
ities for 7 = A and j = B are after measuring each 
particle. These will then serve as the prior probabili- 
ties (J^^^"* and qg'^^'^ for the next measurement, and the 
overall success probability will be the probability of ob- 
taining the correct result at the very last measurement. 
The outcome of this measurement will be used as the 
indicator of what the initial state was. Similar to [23l |. 
we show in Appendix |B] that one of the posterior prob- 
abilities ^7''^^^ will be equal to the success probability 
P, of the fc'th measurement, while the other will be the 
error probability P, = 1 — P,. We then know that ei- 

= 1 — P, if the outcome 



ther g(f+i' 



P, and 



e'^'^ was obtained, or vice versa if the outcome i^j^' 
was obtained. To simplify the notation in the following. 



we define 



Combining 



cos^ 9, = 

Eqs. (|57|) and ([50)1 we can then establish the recursive 
relation 



i?fe = Vl-4Pfc-i(l-Pfc-i)4 



whose solution is 



Rk = 



\ 



(38) 



(39) 



1=1 



From this we see that the probability of obtaining the 
correct result when measuring particle number n, the last 



1 1 



(40) 



fe=i 



where qA = q^"^ and qs = qg'^ are the priors before the 
start of the whole measurement series. When we ap- 
ply this to measuring cat size, we assume equal weight 
for the two branches, so that qA ^ IB = 1/2, and 
Pn = 1/2 + l/2yjl — Y\, cj.. Now if we employ the same 
reasoning as went into deriving Eq. (jlOp for the success 
probability of the optimal collective rt-particle measure- 
ment, we easily obtain that this is identical to P„ in 
Eq. (|40)) . Hence, when the branches are product states, a 
sequence of single-particle measurements with Bayesian 
updating has the same success probability as the optimal 
n-particle measurement. 

The above discussion was carried out entirely in terms 
of distinguishable particles. The result generalizes partly 
to bosonic system, but not entirely. The result holds if 
each of the branches are single-mode Fock states with all 
N particles in the same mode, i.e. \^ a) = (a^)^ |0) /N\ 
and \^ b) = (^^)^ |0) where the modes created by 

and are not necessarily orthogonal. If we then write 

= cos 9c\+ sin 9 c^ and = cos cj, — sin 9 c\^ in anal- 
ogy with Eq. , where cj, and c|, are creation operators 
for orthogonal modes x and y, the bosonic rt-RDMs that 
we obtain using the techniques from Section Hill are iden- 
tical to those we obtain for distinguishable particles using 
Eqs. and ([M)) . Furthermore, the action of the opti- 
mal measurements obtained in the bosonic case can (at 
least in principle) be realized through Kraus operators 
consisting of a single annihilation operator for each mea- 
surement. This simply annihilates a single boson without 
changing the joint state of the system in any other way. 
Hence all conclusions obtained for distinguishable parti- 
cles carry over to the bosonic case in this situation. 

However, if each branch is a more general Fock state 
with more than one occupied mode, i.e. of the form 



\A)^\{al |0) 



B)^\{bl\Q) (41) 



where a\ and a\, may create particles in different modes 
(not necessarily orthogonal) when k ^ k' , then the single- 
particle Bayesian updating measurement protocol de- 
rived above for distinguishable partcles cannot even be 
implemented. Since the particles are not distinguishable 
and cannot be addressed individually, there is no way to 
associate a single value of k with each measurement, and 
hence no way to optimize each single-particle measure- 
ment in the way we did above. Furthermore, if the modes 
associated with different a|, or are not orthogonal, then 



the branches in Eq. (|4T|) in fact contain entanglement be- 
tween modes, and measuring one particle will therefore 
change the state of the remaining system and affect sub- 
sequent measurements. Hence the protocol described in 
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this section only works for bosonic systems when each 
branch is a Fock state with aU particles in a single mode. 



IV. n-RDM ENTROPIES AND RELATED 
CATTINESS MEASURES 

We now analyze the von Neumann entropy of the n- 
RDM and show that this provides insight into how 
meaningful it is to treat the state |^) of Eq. (|13p as a 
two-branch cat state. Calculating the entropy of the n- 
RDM also allows us to compare our cat-size measure to 
an earlier one, the so-called "disconnectivity" introduced 
by Leggett [2g |. 

The von Neumann entropy of a density matrix p is 
given as 

5 = -tr[plnp] = -^panp., (42) 

i 

where {pi} are the eigenvalues of p. Analogous to the 
Shannon entropy of a probability distribution, this quan- 
tity tells us how much information is encoded in the 
knowledge of the physical system represented by the den- 
sity matrix. Equivalently, it can be viewed as the min- 
imum amount of ignorance, we can have about the out- 
come of any measurement on a system represented by a 
given density matrix, where the minimization is over all 
possible measurements encompassed by the density ma- 
trix, i.e., all possible rt-particle measurements in the case 
of an n-RDM. To evaluate the von Neumann entropy Sn 
characterizing n-particle measurements on a cat state we 
need the n-RDM of the fuU state l^*) and not just 
that of the individual branches. This is calculated for the 
states of Eq. (fT3| in Appendix [X] 

Before analyzing the entropy of p^"' for Eq. (fT5|) . we 
first summarize how the entropy should scale for general 
classes of cat-like and non-cat-like states. In general, for 
an experiment that has d equally likely outcomes, the en- 
tropy of the probability distribution is simply In d. If not 
all outcomes are equally likely, then the entropy S will be 
less than Ind. Therefore, if the probability distribution 
of a measurement has entropy S, then the measurement 
must have at least distinct outcomes. This further 
means that, since the von Neumann entropy of a density 
matrix is the minimum entropy of any measurement de- 
scribable by that density matrix, any measurement on a 
system whose von Neumann entropy is S must also have 
at least e"^ distinguishable outcomes. 

For a perfect cat state, schematically of the form 

IV') = l/%/2(|0>^^ + |l)^^) with (0| 1) = 0, this im- 
plies that the n-RDM of the system will have a von Neu- 
mann entropy Sn = ln2, independent of n, until n — N 
where S'tv = 0. If we make a single-particle measurent 
in the {|0) , |1)} basis, the outcomes |0) and |1) are both 
equally likely, so the entropy of that measurement is In 2. 
Unless we measure all N particles however, measuring 
more particles gives us no additional information, since 



measuring just one particle completely collapses the sys- 
tem into one of its branches, and hence the entropy of the 
n-RDM for all n < is equal to In 2. For a "poor" cat 
state, e.g., one of the form |^) ~ I'/'o)'^^ + l*^!)*^^ with 
((^0 \ <t'i) 7^ 0, we cannot distinguish the two branches 
perfectly with an n-particle measurement. One can show 
that the von Neumann entropy in this case will be less 
than In 2. However, as we measure more and more parti- 
cles, the branches become more and more distinguishable 
as they approach orthogonality in the thermodynamic 
limit. Hence the von Neumann entropy will asymptoti- 
cally approach In 2 as n grows. It will then decrease to 
zero again, in a symmetric fashion, as n approaches A^, 
as more and more information about the coherence of the 
branches becomes available. 

Unlike such cat-like states, the entropy of the n-RDM 
of completely generic (pure) states will usually not level 
out as n increases. For a generic state, measuring n par- 
ticles is not likely to tell us very much about the effect of 
adding an n -|~ I'th particle to the measurement. There- 
fore, the number of distinguishable outcomes will usu- 
ally keep increasing with n, until it reaches ^ N/2. At 
that point, we will start gaining enough phase informa- 
tion that the entropy will start decreasing again. At this 
point, the number of particles that we are tracing out be- 
comes smaller than the number of particles we are keep- 
ing, so the entropy can increase no further, and instead 
drops steadily, until it reaches zero at n = A^ (in a pure 
state). 

We turn now to the entropy of p*^"^ for the Gaussian 
cat states defined by Eqs. and (HH). This is plotted 
as a function of n for various parameter combinations 
^0 and a in Figure [31 under the simplifying assumption 
that N ^ n (since we restrict ourselves to this region, 
the drop in entropy as n N cannot be seen). As ex- 
pected from the above general arguments, when a = 
the entropies asymptotically approach In 2 as n ^ oo . 
This means that as we measure more and more particles, 
there exists a von Neumann measurement with exactly 
two distinguishable and equally likely outcomes. In con- 
trast, for (T > the entropy of the n-RDM seems to 
grow without any upper bound, in an approximately log- 
arithmic fashion. This means that, regardless of what 
kind of n-particle von Neumann measurement we make, 
as n — > oo there will always be an ever increasing num- 
ber of distinguishable outcomes. Our state is hence not 
just branching into a nice cat with two cleanly distinct 
branches, but instead developing a whole canopy! This 
canopy keeps growing with n. Hence it does not really 
make sense to view 1 4') as any kind of two- branch or even 
a d-branch cat state in this situation. Instead, it is sim- 
ply some more complicated kind of generic superposition 
state. (The zigzag-pattern for large values of a is caused 
by the factor of (-1)"+'=-' in Eq. (|IIT|) in Appendix O 
which results in a different behaviour for even and odd 
values of n when cr ^ 0, due to interference between con- 
tributions with a given 9 in Eq. (|13p for odd values of 
n.) 
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The von Neumann entropy of the density matrix of the 
fuU state I^P) has previously been used to define a measure 
of cat size referred to as the disconnectivity D by Leggett 
[i,[2i|. To compute D, the entropy Sn of the n-RDM is 
calculated for successively larger n. For each n one also 
finds the minimum total entropy of any partition of the 
n particles, i.e. min,„ (5„, + Sn-m), where the minimum 
is taken over all m from 1 to n — 1. One then defines the 
ratio 113 



mmi<„i 

and the disconnectivity of the system, D, is defined as 
the highest integer n for which /?„ is smaller than some 
"small" fraction <C 1 (/3i is defined to be 0). Thus 
D — max(n|/3„ <C 1). The motivation for this measure 
is that as long as n is smaller than the total number 
of particles needed to observe perfectly the coherence of 
the joint state of all iV particles, the entropy Sn will be 
nonzero since some information about the coherence is 



being neglected when N — n particles are being traced 
out. Subdividing the system further will only neglect 
more information and increase the total entropy, so that 
Sm + Sn-m > Sn and /3„ < 1. As n approaches the num- 
ber of particles sufficient to capture the full coherence of 
the system, Sn and thus /3„ will approach zero. How- 
ever, if n can increase further beyond this point, then 
the denominator will also vanish, and /3„ jumps again 
to 1. Thus the first value of n at which all coherence is 
taken into account will be the largest number for which 
/?„ <C 1. The term "coherence" is used here quite gen- 
erally in the sense of correlations. If the system is made 
up of distinguishable particles and in a pure state, then 
these correlations will be equivalent to entanglement and 
the entropy Sn of the n-RDMs is identical to the bipartite 
entanglement entropy between the n particles included in 
the n-RDM and the N — n particles being traced out [2l[ . 
However, for indistinguishable particles, definition of en- 
tanglement must be made with care, since states with 
little "useful" entanglement can still look very entangled 
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FIG. 4: (Color online.) n-RDM entropies Sn and disconnectivity ratios f5„ for finite systems with A'^ bosons described by 
wavefunctions of the form of Eq. (|13p with three different sets of parameter values 6o and a. The disconnectivity D is defined 
as the largest n value for which /3n <C 1 and is clearly equal to in all three cases shown here. The vertical dashed line 
indicates cat size Co.oi for precision 5 — 0.01. The lower right plot shows the effective cat size as a function of 5 for the three 
cases. 



if one views single particles as good subsystems, due to 
the requirement that the total iV-particle wavefunction 
be symmetrized or anti-symmetrizcd with respect to per- 
mutation of particles [la, [l^, 11^, 22} ■ We will comment 
on these issues in more detail below. 

There are both similarities and important differences 
between the disconnectivity and our measure of effective 
cat size. Both are based on considering how many par- 
ticles must be measured to obtain a specific kind of in- 
formation about the state or its components. But while 
Cs asks how many particles must be measured to dif- 
ferentiate between the two branches composing the to- 
tal state, the disconnectivity D asks how many particles 
must be measured in order to observe all or nearly all 
correlations in the full quantum state. It also does not 
address whether or not the state is naturally divided into 



branches. These differences are reflected in very differ- 
ent numerical results. For the bosonic systems treated 
above, where we have assumed that N is large and made 
approximations based on n 't^ N (see Appendix |^ for a 
full description) , explicit calculation for a range of a and 
6q values shows that Sn increases monotonically with n 
for the whole range treated (except for some minor os- 
cillations between odd and even values of n), so that /3„ 
does not drop below 1/2 until the assumption n N 
is no longer valid. This means that the disconnectivity 
must be of order A'^ for all parameter values a and Oq- 
In contrast. Figure [5] shows that for all values of a our 
measurement-based measure can give values of effective 
cat size Cs much smaller than N, depending on the value 
of 9o. 

In order to make a more direct comparison of D with 
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Ca, we have also calculated Sn for n from 1 to for 
a finite value of N and used this to evaluate the dis- 
connectivity directly for some specific examples. We use 
9o = 97r/40, (T = to study a system close to the full 
overlap situation (^o = 7r/4). We use two examples at 
6*0 = tt/8 {a = 0, N = 100 and a = 7r/16, N = 20) 
for study of an intermediate system and for analysis of 
the effect of nonzero spreading. The n-RDM entropies 
and disconnectivity ratios /3„, Eq. (|43)) . are plotted for 
these three cases in Figure|31 The values of measurement- 
based cat size Cs obtained for these parameters are su- 
perimposed as dashed vertical lines and the bottom right 
panel shows the sensitivity of Cs to the precision S for 
these three cases. It is evident that for all three cases, /3„ 
is more or less constant at a value larger than one half 
and drops to a small fraction substantially smaller than 
this value only at n ^ N. Hence the disconnectivity D 
is equal to or very close to N in all cases. In contrast, 
our cat size measure based on distinguishability gives a 
cat size Cg that is substantially less than TV for all three 
examples. With an error threshold S = 0.01, we obtain 
Co. 01 = N/5 for 6*0 — 7r/8 and cr = or ct = 7r/16, and 
Co. 01 = for Oq ~ 97r/40. Furthermore, the bottom right 
panel shows that in all three cases Cs N for all small 
6, so that our measure differs from disconnectivity for all 
reasonable error thresholds. 

This difference between disconnectivity and 
measurement-based cat size is not totally unexpected. 
In order to observe perfectly all the correlations in the 
states of Eq. [131 one does indeed need to measure all 
or nearly all particles in the system, even when the 
branches are non-orthogonal, unless ^o — 7r/4. However, 
it is clear that except when the branches are orthogonal, 
it is not possible to tell them apart with near certainty 
without measuring more than one particle and one hence 
obtains a reduced effective cat size. Only when we have 
a perfect cat with truly orthogonal branches, e.g., as in 
an ideal GHZ state, will the two measures agree. For 
other states the two measures can be regarded as char- 
acterizing different aspects of the quantum correlations 
in a quantum state. 

Another important aspect of disconnectivity can be 
seen by applying it not to cat-like states but to Fock 
states, i.e. states of the form jvl/) (x a'^ '^fe^^"'^ |0). For 
these states explicit calculation of the n-RDMs and their 
associated entropy Sn shows that D — N for all k ex- 
cept fc = 0, where one obtains D — 1 (see Appendix [C|) . 
In contrast, since Fock states have no branches in the 
second-quantized formalism employed here, n,„in > N 
and the measurement-based cat size measure gives a cat 
size Cs = (see Section [Hi. They also have no entangle- 
ment when expressed in a second-quantized occupation- 
number basis. Thus it may seem puzzling that D can be 
large. However, we note that the disconnectivity relies 
on the entropies of the n-RDMs for its definition, and 
Sn treats individual particles as the fundamental sub- 
systems into which the system is divided and measures 
the correlation between them. As noted in many recent 



papers, this is not appropriate if one is dealing with a 
system of indistinguishable particles, since the system 
can then appear to exhibit full A^-particle entanglement 
simply due to the fact that the wavefunction has to be 
(anti-) symmetrized under exchange of particles. This fic- 
titious entanglement, which has been referred to as "fluffy 
bunny" -entanglement in the literature [33][2§| and which 
goes away if one treats only the modes as good subsys- 
tems instead of particles, is however necessarily present in 
the entropy of the n-RDM, Sn- The fluffy bunny entan- 
glement contribution to disconnectivity is non-zero for all 
states other than those that can be written as Fock states 
with only a single occupied mode. Consequently the dis- 
connectivity of a system of indistinguishable particles will 
be large for all states that are not of this latter special 
kind, whether they are superposition states or not. This 
suggests that one reason for the much larger values of D 
than Cs found here for the states of Eq. (jl3p is inflation 
of the disconnectivity cattiness by fluffy bunny entangle- 
ment. We note that redefining disconnectivity in terms 
of reduced density matrices of modes instead of particles, 
while possible in principle, will however be strongly de- 
pendent on the specific choice of modes. Nevertheless, a 
mode disconnectivity would be limited by the number of 
modes, and for a quantum condensate it is hence likely 
to also be substantially smaller than the total number of 
particles included in the description. 



V. APPLICATION TO CAT STATES OF BEC IN 
A DOUBLE WELL POTENTIAL 

Finally, we apply our measure of cat size to a realistic 
system of bosons in a double well. We consider numerical 
results that have been obtained for bosons with attrac- 
tive interactions in a spherically symmetric 3-dimensional 
harmonic trapping potential, which is split in two by a 
gaussian potential barrier in the xy-plane, forming a dou- 
ble well in the z-direction [3l|. The numerical calcula- 
tions were made using variational path integral Monte 
Carlo (VPI) 32, 33] with 40 interacting bosonic atoms. 
The Hamiltonian used was 



--\7~ + -rf H -=== 

2 ' 2 ' \/ATTa^ 



e 2<T 



+ ^^nt(r,-r,) (44) 

where the sums run over the coordinates of each of 
the 40 atoms and VJ, is a variable barrier height for the 
gaussian potential separating the two wells. Energies 
are given in units of tiuj/2, where lo is the frequency of 
the ground state of the harmoni c trapp ing potential, and 
lengths are given in units of ^Jhjmuj. The two-particle 
interaction potential Vint used here was a Lennard- Jones 
potential 



(45) 
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FIG. 5: Particle number distribution between two wells calculated from Monte Carlo simulation for bosons in a double well 
potential (open symbols) compared with best fit distributions given by Eq. (|13p with a Gaussian spread function f{9) (solid 
lines), hlj = 0.15 for all four cases. 



with Lennard- Jones energy i?Lj and length olj- The 
Lcnnard- Jones potential parameters E^j and a^j de- 
termine the scattering length a It thus provides 
a model potential that allows us to design a compu- 
tationally efficient sampling scheme for a given scat- 
tering length, a (3lj |. Formation of cat states require 
a negative value of a. For a realistic cold atom sys- 
tems with attractive effective interactions such as ^Li 
(a = -14.5 A (H), we find that stable cat states can 
be formed with ~ 1000 atoms in a trap of linear dimen- 
sion aiio = h/mu = 13,000 A, using suitable values of 
Lennard- Jones parameters. 

To compare with our model states in Eq. (|13p . the nu- 
merical data was used to find the probability distribution 
P(n) for finding n of the iV = 40 particles on one side 
of the double-well. This was done for three cases with 
Vt = 10, 15, 20 and with E'lj = 50, olj = .15 in all three 
cases, and for one case with Vb = 120, E'lj = 150, and 
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TABLE I: Best fit of 6o, a, and effective cat sizes Co.oi at 5 = 
0.01 and Cjo-i at 5 = 10~*, for four numerically calculated 
distributions of bosons in a double- well potential, olj = 0.15 
in all cases. 



Olj — -15 (the last choice of extra high potential bar- 
rier and strong attractive interaction was made to get as 
close to a maximal cat state as possible) . We then fit the 
probability distribution Pijia) for the number of parti- 
cles in mode a calculated from the states in Eq. (|13p . to 
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the numerically calculated distributions in each case by 
varying 9o and a to obtain the smallest possible difference 
between the two distributions in the least mean square 
sense. The fitting had a resolution of O.IOtt in 6*0 and 
.OOStt in a. The resulting best fit values for each case are 
shown in Table [H along with the effective cat sizes Co.oi 
for S = 0.01 and C10-4 for 6 ~ 10~^ calculated using the 
states of Eq. with the fitted values of 6*0 and <j (the 
numerical precision in the calculations do not warrant 
smaller values of S). The corresponding fitted number 
distributions are compared to the VPI distributions in 
Figure El showing a very good fit for the cases studied 
here. Note that this does not imply that our states give 
the correct phases between the superposed states, since 
we are only fitting to the number distribution. However, 
given that Eq. (|13p with cr = gives the exact ground 
state in the mean- field limit (l7| . it is reasonable to ex- 
pect that Eq. (|13p constitutes a good approximation to 
the true states. Our comparison with the distributions 
calculated from VPI Monte Carlo supports this expecta- 
tion and also implies that the probability distributions 
(but not necessarily the amplitudes [3l|) can be accu- 
rately described by a two-mode approximation. 

TableUshows that for the lowest barrier height Vb — 10 
we do not really get a cat state at all, since the low bar- 
rier height results in large tunneling, which allows the 
particles to overcome their attractive interactions and 
distribute themselves almost binomially between the two 
wells. The best fit value of 9o (0.227r), is less than one 
a away from the complete-overlap value 7r/4, and the ef- 
fective cat size is correspondingly zero since the branches 
are strongly overlapping. As the barrier height Vb is in- 
creased for a given attraction strength i?Lj , the tunneling 
rate decreases, and it becomes more favorable for all par- 
ticles to sit in one well. However, since the tunneling 
amplitude is still finite, the lowest-energy state is not 
a Fock state but rather a superposition state of nearly 
all particles being in either one well or the other, i.e., 
a cat state. Thus, Cg increases with Vb. In the most 
extreme example here, Vb = 120, E^j = 150, the tunnel- 
ing amplitude is extremely small and the branches have 
negligible overlap, resulting in an ideal cat state Cs = 40 
for N = 40. As expected, we see that Cs does depend 
on the value of the precision 6, becoming smaller as 6 
decreases. We also see that the decrease in cat size is 
greater for larger a values, while for the most "catty" 
case (Vf, — 120 and Slj — 150), where a is practically 
zero, Cs is not affected at all by reducing 6 from 10^^ to 
10-'*. 

We also calculated the disconnectivity D for these 
states and find that D = = 40 in all four cases. This 
may appear initially somewhat surprising, especially for 
the case of — 10 (top left panel in Figure O since in 
that state the branches are almost completely overlap- 
ping, and resemble a binomially distributed state more 
than a cat state. However, even in this case, since the 
distribution is not exactly binomial, there must be some 
entanglement between the particles. Furthermore, all A^ 



particles must be involved in this entanglement since they 
are indistinguishable. As discussed in Section lTVl this co- 
herence between all particles leads to a large value for D, 
even though the state cannot be reasonably called a cat 
state in any way. 

VI. CONCLUSIONS AND FUTURE WORK 

We have presented a measure of the effective size of su- 
perposition states in general quantum systems, i.e., the 
number of effective subsystems that can describe the su- 
perposition, that is based on how well measurements can 
distinguish between the different branches of the state. 
Our measure does in general require one to consider 
coherent multi-particle measurements, although we find 
that for the special class of states considered in 8], a 
procedure using only single-particle measurements can 
be useful. The resulting " cat size" measure is dependent 
on the precision to which the branches are to be distin- 
guished. Application of this measurement-based measure 
to generalized superpositions states of bosons in a two- 
mode system predicts cat sizes much smaller than what 
is predicted from the earlier measure of disconnectivity 
that was proposed in [2^ . Analysis of disconnectivity for 
specific examples showed that for indistinguishable par- 
ticles this quantity is large for a much wider variety of 
states than superposition states, including single-branch 
Fock states, due to the inclusion of particle correlations 
induced by (anti-)symmetrization. 

We expect that the new measure will be useful for 
comparing the effective size of superposition states in 
different kinds of physical systems, including those with 
macroscopic numbers of constituents. We have shown 
that the generalized superposition states studied here 
can be fit to realistic numerical simulations of bosons in 
a 3D double-well trapping potential, and have analyzed 
the cattiness of superposition states of these interacting 
bosons as a function of their interaction strength and 
of the barrier height. Future directions include applying 
our measure to more complicated systems that have been 
realized experimentally, in particular to the experiments 
with superconduc ting loops reported in Q and [1]. In a 
very recent paper [36l | , a different cat size measure was de- 
fined and applied to the three-Josephson junction circuit 
reported in [3] , and the cat size according to that measure 
found to be extremely small (of order 1). It would thus be 
of great interest to evaluate the new measurement-based 
measure of cat size for superpositions of superconducting 
loops. 
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APPENDIX A: CALCULATION OF n-PARTICLE 
REDUCED DENSITY MATRICES 



Inner products between the states (/(J^' (9) y are com- 
puted in the c, d basis using standard methods, giving 



N~l 



(0| (e-*^c - ie^'^d)'^ (e'^'c^ + ie-^^'rft) ^ |0) 
1.N (e'C^-'^') + e-'(^-«')) (0| {e"^' c - le^'df^' (e^'' + ie-^<^'SY |0) 

A^cos {e^e')((i,[^'^\e) 



NiN - l)cos^ {9 - 9') (^<p[^'^^ (9) I 4"^-^^ (0')) 



(Al) 



N\ cos^ {9^9')~mJ—S{ 



and similarly 



0^(0)1^^(0.)^^^, ,^3A^(0. 



(0) 



4") 



9') 

9'): 



(A2) 
(A3) 



The i5-function approximations are valid in the limit of 
large A''. We have assumed that 9 + 9' \s bounded to lie 
between ±7r/2. 
Defining 



\Paf})l - 




.-k ^l^n-l 



[N) 



(A4) 



with a/3 — AA, BB, AB or BA {paa and pbb corre- 
spond to /5^^ and p^g^ as defined in llll A|l . and using the 



action of the operators c, d on the branches 
leads to: 



(9) ) and 



,i-\9) 
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(A6) 



(A7) 



(A8) 



where in the last steps we have made use of the above 
delta function approximation. 

Using the gaussian form in Eq. (PTjl for the amplitude 
spreading function, results in the following analytic forms 
for the rt-RDM matrix elements: 



The traces of the two first matrices are already equal to 
1, so no further normalization is necessary. The trace of 
the matrix defined in Eq. (lAlip is given by 




(AlO) 

The n-RDM of the full state 1 5') , which we use for calcu- 
lating entropies in IIV[ requires also the sum of and 
PiiA which is given by 



k=0 

= 2(e ^ — + (-1) e ^ — j (A12) 



~{n) ~{n) 

Pab + Pba 




(All) so that the final form of the symmetrized ?i-RDM p*-"-*, 
properly normalized, is 
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\Paa + Pbb + Pab + PbAj j 

J / ~(n) , ~(n) , ~(n) , ~(n) \ 



2" [2 + 2 (E_ + 

I 



+ 1 g-2(fc-0^ 



(A13) 



where i?± = exp 



- (00 ± ^/4)' /2a2 



APPENDIX B: DERIVATION OF gV 

7 = A OR B 



(fc+i) 



Pfc FOR 



After a good deal of algebra, using Eqs. 



- ([551) . the 



fact that = Ij and moving factors between the 

two sides of Eq. (jB4[) . both sides can be reduced to 

cos^ 20fc (1 - cos^ 20fc) , (B5) 



To show that one of the prior probabilities (^a^b'' of 
branch A or B before performing the (fc + l)'th mea- 
surement in Section IIII Bl will be equal to the success 
probability Pfe of identifying the correct branch in the 
/c'th measurement, first note that using Bayes' theorem 
and the definition of conditional probabilities, the suc- 
cess probability Pk as given by Eq. ((37|) can equivalently 
be written as 



proving that indeed P{A\E^^'^ 



P(B|i?^'^^ Finally, 
since the measurement on particle fc must give either the 



(fe) 



or E^j^^ , we have P{E^ 



outcome E 
so that Eq. (|Bip reduces to 



P{E 



Pk = P{A\E'^;>)^P{B\Ef). 
which is what we wanted to show. 



(B6) 



Pk = P{A\Ef)P{Ef) + P{B\Ef)P{Ef) (Bl) 

where P{-f\E'^^^) is the posterior probability that the 
state is |^'7) given that the measurement on the fc'th 
particle gave the outcome E^^\ and P{E^^) is the to- 
tal probability that the measurement gives the outcome 
Ej , irrespective of what the state is. Second, note that 
P(A\E^^^) and P{B\E^j^'>) arc in fact equal. This follows 
from 



P{Ef) 



(k) 



q^'^PiE^^'^lB) 
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gt,"'p(£^'°'|J3) 
q(J;->PiE^^-> lA) 



and similarly 



P{Ef\B)qf 
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q^^^P(E^^^\B) 



These will be equal iff 
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qfP{E's'\A) 
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q^^^P{E'^>\B 
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APPENDIX C: DISCONNECTIVITY OF FOCK 
STATES 

In this appendix we show that the disconnectivity, D, 
determined by Eq. is equal to the total particle num- 
ber N for all Fock states that have more than one mode 
with non-zero occupation number. 

A Fock state in a second-quantized system with d 
modes, occupation numbers n = (m, 712, . . . , n^) and a 
total of N particles has the form 



k=l 



|0) 



(CI) 



with rik — N . We assume here that the particles 
are bosons, although this does not affect our final con- 
clusion. We then define a symmetrized n-RDM p'"-* 
by generalizing Eqs. ([^0]) and For this we use 

p = (pi,P2, ■■■,Pd) and q = (gi, 92, • ■ • , gd) as upper and 
lower indices, representing the number of creation and 
annihilation operators, respectively. 



(„)^P_ {N-n)\ 



N\ 



^JUkPkUk 



-.Id 



(C2) 



subject to the constraint that "^kPk ~ l^k Ik ~ n. For 
a Fock state Eq. (jC2p is non-zero only for p = q, i.e. the 
n-RDM is diagonal. Furthermore, we must have Pk,qk 
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rife for a given matrix element not to vanish. For the 
case N = n, the only non-zero matrix element is then 
p = q = n, i.e. the A^-RDM p'^^ has only a single 
matrix element equal to 1 on the diagonal and the rest 
are equal to zero. Hence the entropy is Sn — 0. On 
the other hand, if n < A'^ and if there is more than one 
rife > 0, there will be at least two different p = q for 
which ^ 0, so that p*^"^ must have more than 

one non-zero eigenvalue. Therefore 5„ > for all n < 
N. This implies that the numerator of f3]y in Eq. ([^5]) 



vanishes while the denominator does not. Hence Pn = 0, 
so that n = iV is the largest n for which /3„ ^ 1, and 
consequently the disconnectivity is D ^ N, provided that 
there is more than one mode with non-zero occupation 
number. If only one mode is occupied, p^"-* has only 
a single non-zero eigenvalue (equal to 1) for all n, and 
therefore /?„ = 1 for all n > 1. Since /3i = by definition, 
we therefore have D = 1 for a Fock state in which only 
a single mode is occupied. 
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